#Check the second order code model
@testset "test misocp gf" begin
    @testset "gaslib 582 case" begin
        println("gaslib 582 - MISCOP")                                        
        result = run_gf("../test/data/gaslib-582.json", MISOCPGasModel, misocp_solver)
        @test result["status"] == :LocalOptimal || result["status"] == :Optimal
        @test isapprox(result["objective"], 0; atol = 1e-6)
    end
end


# Check the full nonlinear model
@testset "test minlp gf" begin
        @testset "gaslib 40 case" begin
            println("gaslib 40 - MINLP")                            
            result = run_gf("../test/data/gaslib-40.json", MINLPGasModel, minlp_solver)        
            @test result["status"] == :LocalOptimal || result["status"] == :Optimal
            @test isapprox(result["objective"], 0; atol = 1e-6)
        end
        @testset "gaslib 135 case" begin
            println("gaslib 135 - MINLP")                                    
            result = run_gf("../test/data/gaslib-135.json", MINLPGasModel, minlp_solver)        
            @test result["status"] == :LocalOptimal || result["status"] == :Optimal
            @test isapprox(result["objective"], 0; atol = 1e-6)
        end    
end






